function m = model_moment(x,f,fun)

    m = sum(fun(x).*f)/sum(f);

end